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1 MOTIVATION 



ABSTRACT 

One of the main channels of interactions in galactic nuclei between stars and the central 
massive black hole (MBH) is the gradual inspiral of compact remnants into the MBH 
due to the emission of gravitational radiation. This process is known as an "Extreme 
Mass Ratio Inspiral" (EMRI). Previous works about the estimation of how many 
events space observatories such as LISA will be able to observe during its operational 
time differ in orders of magnitude, due to the complexity of the problem. Nevertheless, 
a common result to all investigations is that the possibility that a compact object 
merges with the MBH after only one intense burst of GWs is much more likely than 
a slow adiabatic inspiral, an EMRI. The later is referred to as a "plunge" because the 
compact object dives into the MBH, crosses the horizon and is lost as a probe of strong 
gravity for eLISA. The event rates for plunges are orders of magnitude larger than slow 
inspirals. On the other hand, nature MBH's are most likely Kerr and the magnitude of 
the spin has been sized up to be high. We calculate the number of periapsis passages 
that a compact object set on to an extremely radial orbit goes through before being 
actually swallowed by the Kerr MBH and we then translate it into an event rate for 
a LIS A- like observatory, such as the proposed ESA mission eLISA/NGO. We prove 
that a "plunging" compact object is conceptually indistinguishable from an adiabatic, 
slow inspiral; plunges spend on average up to hundred of thousands of cycles in the 
bandwidth of the detector for a two years mission. This has an important impact 
on the event rate, enhancing in some cases significantly, depending on the spin of the 
MBH and the inclination: If the orbit of the EMRI is prograde, the effective size of the 
MBH becomes smaller the larger the spin is, whilst if retrograde, it becomes bigger. 
However, this situation is not symmetric, resulting in an effective enhancement of the 
rates. The effect of vectorial resonant relaxation on the sense of the orbit does not 
affect the enhancement. Moreover, it has been recently proved that the production of 
low-eccentricity EMRIs is severely blocked by the presence of a blockade in the rate 
at which orbital angular momenta change takes place. This is the result of relativistic 
precession on to the stellar potential torques and hence affects EMRIs originating 
via resonant relaxation at distances of about ~ 10~ 2 pc from the MBH. Since high- 
eccentricity EMRIs are a result of two-body relaxation, they are not affected by this 
phenomenon. Therefore we predict that eLISA EMRI event rates will be dominated 
by high-eccentricity binaries, as we present here. 



We know, mostly through high-resolution observations of 
the kinematics of stars and gas, that most, if not all, nearby 
bright galaxies harbour a dark, massive, compact object at 
their centres. ( jFerrarese fc Ford|2005{|Kormendy|2004| ). The 
most spectacular case is our own galaxy, the Milky Way. By 
tracking and interpreting the stellar dynamics at the centre 
of our galaxy, we have the most well-established evidence for 
the existence of a massive black hole (MBH) (Eisenhauer 
et al.||2005| iGhez et al.||2005| 120081 iGillessen et al.||2009|) 



Observations of other galaxies indicate that the masses of 
MBH can reach a few billion solar masses (M ). The exis- 
tence of such a MBH population in the present-day universe 
is strongly supported by Soltan's argument that the aver- 
age mass density of these MBHs agrees with expectations 
from integrated luminosity of quasars (Soltan 1982, Yu fc| 
Tremaine|2002[ ) . Many correlations linking the MBH's mass 
and overall properties of its host spheroid (bulge or ellipti- 
cal galaxy) have been discovered. The tightest are with the 
spheroid mass (I Haring fc Rix||2004 l, its velocity dispersion 
(M — a relation, Tremaine et al.|2002 I and degree of concen- 
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tration (E rwin et al.|2004[). Consequently, understanding the an orbit for which t GW <C (1 — e)i r (|Amaro-Seoane et al 



origin and evolution of these MBHs necessitates their study 
in the context of their surrounding stellar systems. 

The ideal probe of these regions is the gravitational ra- 
diation (GWs) that is emitted by some compact stars very 
close to the black holes, and which will be surveyed by 
cLISA/NGO (evolved Laser Interferometer Space Antenna / 
New Gravitational Wave Observatory [Amaro-Scoanc et ah] 
2012a|b| . This mission will scrutinise the range of masses 



fundamental to the understanding of the origin and growth 
of supermassive black holes; i.e. MBHs with masses below 
10 7 M o . 



2 EMRIS AND DIRECT PLUNGES 

For a binary of a MBH and a stellar black hole to be in 
a LISA-like band, it has to have a frequency of between 
roughly 1 and 10 Hz. The emission of GWs is more ef- 
ficient as they approach the LSO, so that the observatory 
will detect the sources when they are close to the LSO line. 
The total mass required to observe systems with frequen- 
cies between 0.1 Hz and 10~ 4 Hz is of 10 4 - 1O 7 M . For 
masses larger than 10 7 M Q the frequencies close to the LSO 
will be too low, so that their detection will be very difficult. 
On the other hand, for a total mass of less than 10 3 M Q we 
could in principal detect them at an early stage, but then 
the amplitude of the GW would be rather low. 

To interact closely with the central MBH, stars have to 
find themselves on "loss-cone" orbits, which are orbits elon- 



gated enough to have a very close-in periapsis (Frank & Rees 



1976 



Lightman & Shapiro 1977; Amaro-Seoane & Spurzem 



2001). The rate of tidal disruptions can be established semi- 



/ analytically if the phase space distribution of stars around 



the MBH is known (Magorrian & Tremaine 1999 Syer & 
|Ulmer||l99"9"l |Wang fc Merritt||2004[ for estimates in models 
of observed nearby nuclei). To account for the complex in- 
fluence of mass segregation, collisions and the evolution of 
the nucleus over billions of years, detailed numerical simu- 
lations are required, however ( David et al.|1987a|b Murphy 



et al.||1991| |Freitag fc Benz||2002l |Baumgardt et al.| 
Amaro-Seoane|2004| |Freitag et al.|2006||Khalisi et al 



2004 



2007 



Preto fc Amaro-Seoane|2010||Amaro-Seoane fc Preto|201l[ ). 



As the star spirals down towards the MBH, it has many 
opportunities to be deflected back by two-body encounters 
onto a "safer orbit" (Alexander & Hopman 2003 Amaro- 



Seoane et al.|2007 l, hence even the definition of a loss cone 
is not straightforward. Once again, the problem is com- 
pounded by the effects of mass segregation and resonant 
relaxation, to mention two main complications. As a re- 
sult, considerable uncertainties are attached to the (semi- 
analytical predictions of capture rates and orbital parame- 
ters of EMRIs. 

Naively one could assume that the inspiral time is domi- 
nated by GW emission and that if this is shorter than a Hub- 
ble time, the compact object will become an EMRI. This is 
wrong, because one has to take into account the relaxation 
of the stellar system. Whilst it certainly can increase the 
eccentricity of the compact object, it can also perturb the 
orbit and circularise it, so that the required time to inspiral 
in, t GW , becomes larger than a Hubble time. The condition 
for the small compact object to be an EMRI is that it is on 



2007), with t r the local relaxation time. When the binary 



has a semi-major axis for which the condition is not ful- 
filled, the small compact object will have to be already on a 
so-called "plunging orbit", with e ^ e p i UIlgc = 1 — 4_R Schw /a, 
where -Rgchw is f ne Schwarzschild radius of the MBH, i.e. 
7? Schw = 2GM./c 2 , with M. being the MBH mass. It has 
been claimed a number of times by different authors that 
this would result in a too short burst of gravitational radi- 
ation which could only be detected if it was originated in 
our own Galactic Center ( Hopman et al.|2007 1, because one 
needs a coherent integration of some few thousands repeated 
passages through the periapsis in the eLISA bandwidth. 

Therefore, such "plunging" objects would then be lost 
for the GW signal, since they would be plunging "directly" 
through the horizon of the MBH and only a final burst of 
GWs would be emitted, and such burst would be (i) very 
difficult to recover, since the very short signal would be 
buried in a sea of instrumental and confusion noise and (ii) 
the information contained in the signal would be practically 
nil. There has been some work on the detectability of such 
bursts jRubbo et al.|2006||Hopman et al.|2007||Yunes et~aL] 



2008 Berry fc Gair 2012 1 , but they would only be detectable 



in our galaxy or in the close neighborhood, but the event 
rates are rather low, even in the most optimistic scenarios. 

The typical size of the central MBH can be associated 
with the gravitational radius (radial horizon location) of the 
MBH, which in the case of the Milky Way MBH corresponds 
to approximately i?g chw ~ 1.3 x 10 7 km ~ 4.1 x 10 -7 pc 
(neglecting spin contributions). This number gives a good 
indication of how small are these MBHs in size, which means 
they have a small cross section and hence, the chances a star 
has to "plunge" through the horizon of the MBH directly are 
very small. 

To quantify the probability for star absorption by a 
MBH it is crucial to take into account the location of the 
Last Stable Orbit (LSO) of a test massive body in terms 
of the MBH parameters. According to General Relativity, 
these are the mass M, and its intrinsic angular momen- 
tum or spin S, = a, M.c, where a, is the spin parame- 
ter with length dimension and subject to the constraints 
«: a.c 2 /(GM.) sC 1. The LSO location is given by a sur- 
face in the orbital configuration space that can be described 
in terms of the parameters (p, e, i), where p is the dimen- 
sionless semilatus rectum, e is the eccentricity, and i is the 
orbital inclination (with respect to the MBH spin axis) of 
the orbit. Equivalently one can also use the semimajor axis 
a, or the periapsis location r pcri , instead of p. The periapsis 
and apoapsis radii are given then by: 



GM.p 

z 2 (l + e) 
GM.p 

9(1 - e) 



(1) 



It is well-known (see e.g. Bardeen|1970 1 that the LSO, 
for the case of circular orbits, lies at 3 -Rschw f° r non-spinning 
MBHs, while it is shifted out to 9GM,/c 2 for retrograde or- 
bits and down to GM./c 2 for prograde orbits, and these 
values correspond to the case of "extremal" MBHs, charac- 
terized by maximal spins, i.e. a,c 2 / (GM,) = 1 . Despite this 
fact, traditional EMRI event rate estimations are based on 
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considerations that neglect the spin of the MBH. Taking into 
account the spin one would expect (considering the asym- 
metry between prograde and retrograde orbits) an increase 
of the number of EMRIs since some of the traditionally ne- 
glected "direct plunges" can actually be disguised EMRIs. 



3 ORBITAL GEODESIC MOTION AROUND A 
KERR MBH 

In order to show the importance of the effect of the spin in 
the estimation of the number of EMRIs that will produce a 
significant amount of GW cycles in the band of eLISA, we 
present here two types of calculations. The first one is to 
adapt known results about the stability of orbits of massive 
objects around a MBH to our discussion. The second is an 
estimation of the number of cycles for for orbits which would 
be plunging orbits for a Schwarzschild MBH, or orbits with 
no sufficient cycles when the MBH was assumed to be non- 
spinning for the case with spin. We show that a significant 
fraction of them are actually EMRIs with sufficient cycles to 
be detected by a space-based observatory like eLISA. Parts 
of these calculations, mainly due to the high eccentricities 
involved, require numerical computations. 

At this point it is useful to review some basic character- 
istics of the orbital geodesic motion of massive bodies around 
a Kerr MBH. First of all, the geometry of a Kerr MBH 
is axisymmetric (with respect to the spin axis) instead of 
spherically-symmetric as in the case of Schwarzschild MBHs, 
and this means that the inclination of the orbit with respect 
to the spin axis, t, plays an important role in the dynamics. 
Actually, orbits outside the equatorial plane are not planar, 
like in Keplerian motion or orbits around a Schwarzschild 
MBH, but instead they would precess around a plane with 
a certain inclination i with a frequency that we call fg, 
where 6 refers to the polar Boyer-Lindquist coordinate of 
the MBH ( |Boyer fc Lindquist|p67| |Misner et!j^[T973l ) . In 
addition, relativistic effects cause precession of the periap- 
sis, and this already happens for Schwarzschild MBHs, so 
that we have to consider two more frequencies, f r and f v . 
The first one, f r , is associated with the radial motion and 
the time to go from periapsis, r. 



stant per unit mass square, C, which is associated with an 
extra symmetry of the Kerr geometry, similar to what hap- 



pens in certain axisymmetric Newtonian potentials (Binncy 
fc Tremaine|1987 l. Actually, the set of constants (E, L 2 , C) 



parametrizes the geodesic orbit in the same way as the set of 
orbital parameters (p, e, i) does. Therefore, there is a map- 
ping between these two sets (see Schmidt 2002 I which is 



going to be crucial in the calculations that we present here. 
The explicit form of this mapping is quite complex and we 
do not include it here, we just mention that we used the 



implementation described in Sopuerta & Yunes (20111. 

In the case of a non-spinning Schwarzschild MBH, 
where i and C do not play any role, the mapping is much 
more simple and is given by: 



E 2 _{p- 2 - 2e)(p- 2 + 2e) 



(2) 



(3) 



c 2 p (p — 3 

L 2 G 2 M 2 p 2 
z c 2 (p — 3 — e 2 ) 

Using the symmetries of the geometry of a Kerr MBH 
we can separate the equations for geodesic orbital motion 
so that the trajectory of a massive body, described in terms 
of Boyer-Lindquist coordinates {t, r, 9, <p}, can be written as 
follows 



c c 



r dr A 

dr x ~ 
rf7 



. 2\ E L z 
+ a. ) a. — 

' c c 



+ r A = R(r) 



C L z in 2 i . 

-3 - ^ cot 9 ~ a > I 1 



E^ 2t 
cos ( 



dtp 1 
d7 = A 



E L, A — a.sm 
2a,r m —r+ — 

c c 



Sill 



2 
2 

2 



(4) 

(5) 
(6) 
(7) 



to apoapsis, r apo , and 



In the last set of equations we have introduced the following 
definitions: 

_ GM. 



motion around the spin axis and the time to complete a full 
turn (2 7r) around this axis, or in other words, the time for 
the azimuthal angle tp to increase 2n radians. 

In summary, generic bound motion around a Kerr MBH 
exhibits three fundamental frequencies, (/,., f s , f v ) and this 
implies that the GW emission of an EMRI will be quite 
rich in structure (not only these GWs will contain features 
with these frequencies but also with a number of harmon- 
ics of them), encoding the detailed geometry of the central 
Kerr MBH. The GW emission will back react on to the sys- 
tem and this translates in particular in a change of the or- 
bital parameters (p, e, t) of the orbit. These changes can 
be estimated by considering the energy and angular mo- 
mentum carried away from the extreme-mass-ratio binary 
by the GWs emitted. More specifically, the GW emission 
changes the constants of motion of the geodesic motion, 
namely the energy per unit mass (normalized with respect 
to the star mass, m), E, the angular momentum along the 
spin axis per unit mass, L z , and the so-called Carter con- 



Q = C+(L z -a.Ey , 

2 2 . 2 2/i 

p =r + a. cos o 

A = r 2 - 2r.r + a 2 , = r 2 f + a 2 ., with f = l- — . (8) 

For convenience, we also define the quantity E 2 = (r 2 + 
a\) 2 — a^A sin 2 8 . The first equation tells us how the co- 
ordinate time t changes with respect to the proper time r 
and the other three describe the trajectory in space. One 
can combine the four equations to obtain the spatial trajec- 
tory in terms of coordinate time t, the time of observers at 
infinity, i.e. (r(t), 0(t), <p(t)). 



4 KERR AND SCHWARZSCHILD 
SEPARATRICES 

In figures[l]and[2]we show plots of the location of the LSO in 
the plane a (pc) - (1 — e), including the Schwarzschild sepa- 
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ratrix between stable and unstable orbits, p — 6 — 2e = (Q for 
both prograde and retrograde orbits and for different values 
of the inclination t. Each plot corresponds to a different value 
of the spin, showing how increasing the spin makes a differ- 
ence in shifting the location of the separatrix between stable 
and unstable orbits, pushing prograde orbits near GM./c 2 
while retrograde orbits are pushed out towards 9GM./c 2 . 
The procedure we have used to build these plots is a stan- 
dard one. Briefly, given a value of the dimensionless spin 
parameter s = a,c 2 /(GM.) and a value of the eccentricity 
e and inclination angle i we apply the following algorithm: 

(i) We start from an initial value of the semilatus rectum 
p which, together with the value of the eccentricity fixes the 
apoapsis and periapsis radii through the equations in 
On the other hand, the inclination of the orbital plane can 
also be described in terms of the polar angle 8, so that (9 inc 
represents the inclination angle and is closely related to l 
(see e.g., |Drasco fc~ Hughcs 2 004 \ . The advantage of using 
6 inc is that it is related in a simple way to the extrema of 9, 
# min satisfying -see equation Q 



(9) 



by inc = sign(L z 



sign(L 2 



with 



+ 1 for prograde orbits , 
— 1 for retrograde orbits . 



(10) 



Then, from the condition of extrema of 9 min and its relation 
to # inc , we can find the value of the Carter constant C in 
terms of the energy E and angular momentum along the 
spin axis L z . 

(ii) From the conditions of extrema of periapsis and 
apoapsis, 



= 0, 



(11) 



[see equation jBJ], and using the expression of C in terms 
of (E,L Z ) from the previous point, we find the values of 
(E,L Z ) (and hence of C too). 

(iii) The radial motion for geodesic orbits around Kerr 
has four extrema, the periapsis and apoapsis locations and 
two more radii, r 3 and r 4 , such that 

r ap o > r pcli ^ r 3 > r 4 . (12) 
Actually, r4 always lies inside the horizon radius, 



(13) 



i.e. r 4 < r H . For any stable orbit, it is obvious that the radial 
motion happens inside the interval [r peri , r apo ]. However, for 
orbital configurations with r peri = r 3 , the potential of the 
MBH changes its shape and the orbits become unstable, 
marking the location of the LSO. In this case, we have that 



1 The relation between p and a is a = (GM,/c 2 )(p/ (1 — e 2 )) . 



dR(r) 
dr 



(14) 



where R(r) denotes the right-hand side of the evolution 
equation for the radial position r (see equation IH|. 

The calculations in this algorithm are done numerically, so 
we check whether this condition is satisfied to some tolerance 
level. In the case it is not satisfied, we use this information to 
prescribe the next value of p and come back to the first point 
in the algorithm. The process is repeated until we identify 
the LSO with the desired accuracy. 



5 NUMBER OF CYCLES 

The second type of relativistic computations that we have 
performed concerns the estimation of the number of cycles 
that certain EMRI orbital configurations that were thought 
to be plunging orbits (or orbits with no sufficient cycles) in 
the case of non-spinning MBHs can spend in a frequency 
regime of / £ [10~ 4 , 1] Hz during their last year(s) of in- 
spiral before plunging into the MBH. This is important to 
assess how many of these EMRIs will have sufficient Signal- 
to-Noise Ratio (SNR) to be detectable. The way in which 
these estimations have been done is the following. We start 
with a certain orbital configuration characterized by the or- 
bital parameters (p( )> e (o)i '(()))• Equivalently, we can char- 
acterise the initial orbital configuration by the constant of 
motions (E^, L z ( \, Cfo))- Hence, the idea is to track the in- 
spiral without actually integrating the equations of geodesic 
motion of section[3]or any other type of equations that follow 
the trajectory. Instead, we picture the inspiral as a sequence 
of geodesic orbits, each of them characterized by orbital pa- 
rameters (p(j), era, ( or equivalently, constants of motion 



C w )) withi = 0, 



N , 

i A * plunge ) 



being iVpi unge the 



final plunging configuration. The transition between each 
orbital configuration is governed by the GW emission. Our 
particular algorithm to follow the inspiral goes as follows 
(our implementation uses the formulae in the appendices 
of Sopuerta & Yuncs 2011 and the formulae derived by |Gair| 
fc Glampedakis||2006[ ): ~ 



(i) Given an orbital configuration (p/f, , e./ j) , t(j) ) and its 
associated constants of motion {E^,L Z ^,C^), we com- 
pute the averaged evolution of the constants of motion, 
(E,L Z ,C), using the formulae of (Gair & Glampedakis 



|2006p , which combine post-Newtonian calculations at the 
2PN order with fits to results for the GW emission based 



on the Teukolsky formalism ( Teukolsky 1973 1 (for details 
see|Hughes l (|2000|l2001[);|Glampedakis et al.| ( [2002"l ); |Drasco 
|fc Hughes||2006| ); [Hughes et al.| < |2005| l). 

(ii) For the given orbital parameters we estimate the ra- 
dial period T r , that is the time to go from apoapsis to periap- 
sis and back to apoapsis (see |Schmidt|2002|[Fujita fc Hikida| 
[2009] for details of this computation. In |Fujita fc Hikida|20W 



there is a typo in one of the relevant formulae for our compu- 
tations fixed in the appendices of Sopuerta fc Yunes|20lT l. 

(iii) We compute the change in the constants of motion. 
To that end, due to the fact that the GW emission by 
an EMRI is relatively weak, we do not consider in general 
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Figure 1. LSO for a MBH of mass 4 X 10 4 M Q and a SBH of mass m. = 10 M for a Kerr MBH of spin s = 0.4 (left) and s = 0.7 
(right). The Schwarzschild separatrix is given as a solid black line. Curves above it correspond to retrograde orbits and inclinations of 
i = 5.72, 22.91, 40.10, 57.29, 74.48 and 89.95° starting from the last value (89.95°). In the left panel we can barely see any difference 
from the different inclinations due to the low value of the spin. 



just one passage through periapsis, but several of them, say 
N?™ 1 . Thus, the change in the constants of motion is: 



AL z (j), AC(j)) — (-E(i), L x (f), C(i)) 



x Nf: 11 x T r 
W r 



(15) 



And from here, 



(-E(i+l))iz(!+l)iC(i+l)) — 0®(i)> -k*(<)> + 

(A£ (l) ,AL 2(l) ,AC (l) ). (16) 

(iv) From these new constants of motion and using 
similar techniques to the ones described above for the 
LSO computation we can find the new orbital parameters 
t( i+1 j). We then go back to the first step and 
follow the inspiral until we detect that the new configura- 
tion corresponds to an unstable orbit corresponding to the 
plunge. We estimate the number of cycles as the number of 
times that the stellar object turned around the spin axis, 

cycles 

In Table [TJ we show some results for a series of inspirals 
whose initial orbital parameters are such that the pair (e,p) 
lies on the separatrix between stable and unstable orbits in 
the case of non-spinning MBHs (i.e. p — 6 + 2e). There- 
fore, these are inspirals that under the assumption that spin 
can be neglected would have no cycles in the eLISA band. 
However, we can see from the values for the number of cy- 
cles in the table that many of those systems actually per- 
form a significant number of cycles, more than sufficient to 
be detectable with good SNR. The number of cycles has 
been associated with N v (the number of times that the az- 



imuthal angle ip advances 2ir) which is usual for binary sys- 
tems. However, as we have discussed above the structure 
of the waveforms from EMRIs is quite rich since they con- 
tain harmonics of three different frequencies. Therefore the 
waveforms have cycles associated with the three frequencies 
(fr, fg, ftp) which makes them quite complex and in principle 
this is good for detectability (assuming we have the correct 
waveform templates). Another fact that is worth mention- 
ing is that the cycles quoted in Table [TJ happen just before 
plunge and take place in the strong field region very near the 
MBH horizon. Then, these cycles should contribute more to 
the SNR than cycles taking place farther away from the 
MBH horizon. Regarding the accuracy of these estimations, 
the main sources of error are the approximations made for 
the radiation-reaction effects, which are based on a post- 
Newtonian expansions and fits to results from black hole per- 
turbation theory. Corrections from higher-order terms will 
introduce corrections to these results that depend on the 
EMRI configuration, but those corrections should not affect 
the magnitude of these numbers. In the Table we quote the 
integers that are closer to the numerical result. 



6 IMPACT ON EVENT RATES 

Only a certain fraction of stars in phase space will come 
close enough to interact with the MBH. These stars are said 
to belong to the "loss-cone" (see e.g. Frank & Rees 1976 
Amaro-Seoane fc Spurzem||2001 \ 



For radii larger than 0.01 pc the main leading mecha- 



nism for producing EMRIs is two-body relaxation ( |Hopman 
& Alexander 2005; Amaro-Seoane ct al. 2007; Amaro-Scoanc 
20121, and this is the region of phase-space in which our 
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analysis is applied with priority, since for a Schwarzschild 
MBH one has just direct plunges. For radii below 0.01 pc 
we note that the enhancement in the EMRI event rate due 
to resonant relaxation predicted by |Hopman fc Alexander| 
( 2005 1 is severely affected by the presence of a blockade in 



the rate at which orbital angular momenta change takes 
place. This so-called "Schwarzschild barrier" is a result of 
the impact of relativistic precession on to the stellar poten- 



tial torques, as recently shown by Merritt et al. ( 2011 1 with a 
few direct-summation N— body simulations expanded with 
a statistical Monte-Carlo study. Indeed, this "Schwarzschild 
barrier" has been corroborated in an independent work with 
a statistical study based on a sample of some 2,500 direct- 



summation N— body simulations by Brem et al. (20121 in 
which the authors include post-Newtonian corrections and 
also, for the first time, the implementation of a solver of 
geodesic equations in the same code. This barrier poses a 
real problem for the production of low-eccentricity EMRIs. 
However, high-eccentricity EMRIs, which had been classified 
until now wrongly of "plunges" do not have this problem, 
since they are a product of pure two-body relaxation. This 
is why we will only focus on two-body relaxation for the 
estimation of the rates. 

The event rate can be hence approximately calculated 



that the density profile follows p oc r~ 7 within the region 
where the gravity of the MBH dominates the gravity of the 
stars, with 7 ranging between 1.75 and 2 for the heavy stel- 



lar components 1 


Peebles 


1972 


Bahcall & Wolf| 


1976 1977 


Amaro-Seoane et al. 2004| Preto et al. 2004| Alexander & 


Hopman|2009; Preto & Amaro-Seoane|2010 Amaro-Seoane 


& Preto||2011|) and see Gurevich ( 


1964 


1 for an interesting 



stars within a radius r is 



n(r) = 



(3- 7 ) (M. 



4 7r \m 

Hence, the number of SBHs within a is 



N.(a)~N 



3-7 



Hence, we have that 



2-7 



dN.(a) = (3-7)^(^1 'lo- 



(19) 



(20) 



A^EMRI 



/' 



dN.(a) 



In (6 



Lc) 



t r (a) 



(17) 



with #lc the loss-cone angle, N,(a) the number of stellar 
black holes (SBHs) within a semi-major axis a and a EMRI 
a maximum radius within which we estimate the event rate 



We assume that the SBHs distribute around the 



central MBH following a power-law cusp of exponent 7, 



6.1 The Schwarzschild case 



Livio 2001 Amaro- 



We know that (see e.g. [Alexander 

Seoane fc Spurzem||2001[ ) ~ Jmax/^LC- Since the loss- 
cone angular momentum can be approximated as J L c 
4GM,/c and J max = y/GM,a, we have that 



2 The authors obtained a similar solution for how electrons dis- 
tribute around a positively charged Coulomb centre. 
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MBH mass 


MBH Spin 


Semimajor axis 


Eccentricity 


Inclination 


Time to plunge 


Time in band 


Cycles in band 


M. (M ) 


s 


a (0 ) (P c ) 


— e (o) 


t (0) (rad) 


T , (yrs) 

plunge \ J ' 


T band (y fS ) 


N v 


5.0 • 10 4 


0.30 


9.57- 10~ 3 


10" 6 


0.00 


3.92 • 10 2 


0.06 


914 


5.0 ■ 10 4 


0.30 


9.57- 10~ 3 


10" 6 


0.70 


3.92 ■ 10 2 


0.06 


625 


5.0 ■ 10 4 


0.90 


9.57- 10" 3 


10" 6 


0.00 


3.92 ■ 10 2 


0.11 


4174 


5.0 • 10 4 


0.90 


9.57- 10" 3 


10" 6 


1.00 


3.92 • 10 2 


0.08 


2646 


1.0 • 10 5 


0.70 


1.91 • 10~ 2 


10" 6 


0.00 


7.87- 10 2 


0.29 


6968 


1.0 ■ 10 5 


0.70 


1.91 ■ 10~ 2 


lO" 6 


1.00 


7.86 ■ 10 2 


0.23 


3411 


1.0 ■ 10 5 


0.99 


1.91 ■ 10" 2 


10" 6 


0.00 


4.70- 10 3 


0.38 


8938 


1.0 • 10 5 


0.99 


1.91 • 10" 2 


10" 6 


0.70 


3.92 • 10 3 


0.32 


7892 


5.0 • 10 5 


0.30 


9.57- 10~ 2 


10" 6 


0.00 


4.31 ■ 10 4 


1.98 


8246 


5.0 ■ 10 5 


0.30 


9.57- 10~ 2 


10" 6 


1.00 


4.31 ■ 10 4 


1.61 


3061 


5.0 • 10 5 


0.95 


9.57- 10" 2 


10" 6 


0.00 


5.10 • 10 4 


2.00 


40093 


5.0 ■ 10 5 


0.95 


9.57- 10~ 2 


10" 6 


1.00 


4.31 ■ 10 4 


2.00 


27600 


1.0 • 10 6 


0.30 


1.91 • 10~ 2 


10~ 5 


0.00 


1.27- 10 4 


1.98 


10943 


1.0 • 10 6 


0.30 


1.91 • 10" 2 


10" 5 


1.00 


1.19 • 10 4 


1.91 


3552 


1.0 ■ 10 6 


0.70 


1.91 ■ 10" 2 


10" 5 


0.00 


1.35 ■ 10 4 


1.99 


51308 


1.0 ■ 10 6 


0.70 


1.91 ■ 10" 2 


10" 5 


1.00 


1.20 ■ 10 4 


1.99 


23291 


1.0 • 10 6 


0.90 


1.91 • 10~ 2 


10" 5 


0.00 


1.40- 10 4 


1.99 


58841 


1.0 • 10 6 


0.90 


1.91 • 10~ 2 


10~ 5 


1.00 


1.17- 10 4 


2.00 


38245 


1.0 ■ 10 6 


0.99 


1.91 • 10" 2 


10~ 5 


0.00 


1.43 • 10 4 


2.00 


61726 


1.0 ■ 10 6 


0.99 


1.91 ■ 10" 2 


10~ 5 


1.00 


1.17- 10 4 


2.00 


47678 


5.0 • 10 6 


0.30 


9.57- 10~ 2 


10" 5 


0.00 


1.44- 10 5 


1.93 


5258 


5.0 • 10 6 


0.30 


9.57- 10~ 2 


10~ 5 


1.00 


1.36 • 10 5 


0.00 





5.0 ■ 10 6 


0.70 


9.57- 10" 2 


10~ 5 


0.00 


1.55 ■ 10 5 


2.00 


40687 


5.0 ■ 10 6 


0.70 


9.57- 10" 2 


10~ 5 


1.00 


1.36 ■ 10 5 


2.00 


14936 


5.0 • 10 6 


0.90 


9.57- 10~ 2 


10" 5 


0.00 


1.61 • 10 5 


2.00 


41369 


5.0 • 10 6 


0.90 


9.57- 10~ 2 


10~ 5 


1.00 


1.35 • 10 5 


1.99 


30695 


1.0 • 10 7 


0.30 


1.91 


10" 6 


0.00 


4.02 • 10 6 


1.99 


3089 


1.0 ■ 10 7 


0.30 


1.91 


10" 6 


1.00 


3.78 • 10 6 


0.00 





1.0 • 10 7 


0.70 


1.91 


10" 6 


0.00 


4.27- 10 e 


2.00 


23425 


1.0 • 10 7 


0.70 


1.91 


10~ 6 


1.00 


3.79 • 10 e 


1.98 


8747 


1.0 ■ 10 7 


0.99 


1.91 ■ 10" 1 


10~ 5 


0.00 


1.44- 10 6 


1.98 


22455 


1.0 • 10 7 


0.99 


1.91 ■ 10" 1 


10~ 5 


1.00 


1.18 • 10 6 


1.99 


28589 


5.0 • 10 7 


0.30 


9.57- 10" 1 


10" 5 


0.00 


1.44- 10 7 


0.00 





5.0 • 10 7 


0.30 


9.57- lO" 1 


10~ 5 


1.00 


1.36 • 10 7 


0.00 





5.0 • 10 7 


0.70 


9.57- 10- 1 


10~ 5 


0.00 


1.55 ■ 10 7 


1.72 


4247 


5.0 • 10 7 


0.70 


9.57- 10" 1 


10~ 5 


1.00 


1.35 ■ 10 7 


0.00 





5.0 • 10 7 


0.99 


9.57- 10" 1 


10" 5 


0.00 


1.65 • 10 7 


1.88 


4422 


5.0 • 10 7 


0.99 


9.57- 10" 1 


10~ 5 


1.00 


1.35 • 10 7 


1.52 


4625 



Table 1. This table shows the main properties of some (prograde) inspirals that initially lie on the separatrix (LSO) of non-spinning 
MBHs and hence they would not be detectable in the eLISA band. The numbers in the first five columns have been already introduced 
in the text. The sixth column gives the time it takes for each inspiral to get to plunge. The seventh column shows how much time it 
spends in band assuming the plunge occurs at the end of the eLISA mission time (assumed to be 2 yrs here). The last column show the 
number of orbital cycles in band (during T band ), defined as the number of times that the azimuthal angle <p advances 2tt during the last 
two years before plunge. The number of GW cycles can be then defined as twice this number. 



8Rsc 



(21) 



We assume also that relaxation is dominated by a single 
stellar black hole (SBH) population, since because of mass 
segregation the most massive objects sink down to the centre 
and the light stars are pushed out from the centre. The 
relaxation time at a distance a is 



t r (a) = t 



7-3/2 



(22) 



with 



t = 0.3389 : 



In A G 2 m 2 n 
Since (see e.g. Amaro-Seoane et al.|2007 l 



n = ^NoRl 
1 GM. 



2 

c 



1+7 



we have that equation ( 23 ) becomes 



t ~ 4.26 



^Rl{GM.y 



(3- 7 )(l + 7) 3 / 2 



In A AT, 



o 



M. 



(23) 

(24) 
(25) 

(26) 
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We now define a critical radius at which the two regimes that 
lead the evolution of the EMRI decouple. The first evolution 
is dominated by relaxational processes, via exchange of E 
and J between SBHs on a capture orbit with stars from the 
surrounding stellar system, while in the second regime the 
evolution of the EMRI is totally dominated by the emission 
of gravitational waves. This is given in figure 1 of |Amar(>] 
Seoane et al. (20071 with their red curve. In other words, 



the line gives us the radius as a function of a at which the 
relaxational time at periapsis is approximately equal to the 



timescale defined by the approximation of Peters ( 1964 1 
Hence, we have to solve 



MBH. We then have to either shrink or enlarge it by a cer- 
tain factor function of the inclination and spin, W(t, s), so 
that the effective pericentre of LSO is 



(l-e)o = W(t, s) 



\GM. 



(32) 



This quantity can be derived from the separatrices of the 
figures in section [3] If the orbit can get closer to the MBH 
in the Kerr case, then W(t, s) < 1; otherwise VV(t, s) > 
1. Since the separatrices are nearly parallel, we hence can 
define W(l, s) like 



t r (a)(l - e) = Kt GW (a, e) 
J(a, e) = Jlc 

v 8GM, 

0--e)a= — -= — 



(27) 



In the first equality, K is a factor of order unity. In the 
last equality we assume a Schwarzschild radius and we as- 
sume that the LSO is at 4 x - Rschw Approximating e ~ 1, 
the function /(e) from |Peters| |l964| /(e) = 425/(768^2). 
Hence 



, , , „ /k 24 c 5 a 4 7/2 



(28) 



And so, finally from equations (271, (221 and (26 1 and solv- 
ing for Hemrii we have that 



a EMRI — Rq 



16.97 K (3 - 7) (1 + 7 ) 3/2 In A N, -j^- 



(29) 



Or, absorbing some constants into a newly defined 



W(t, s) 



a LSO, Kerr 
ffl LSO, Schw 



N ^ 



l LSO 



Kerr v^i) 



LSO. Schw 



(33) 



That is, given a spin and an inclination, we take for N val- 
ues of the eccentricity the semimajor axis corresponding to 
the Kerr value, Olso Kerr( e ) an d to the Schwarzschild case, 
a LSO Schw an d we sum for the ratio of the two semi-majors 
over all eccentricities. This allows us to calculate by how 
much the LSO has been "shifted" . 

If w e re do the derivation of (31 1 taking into account 
equation 32 we can obtain AT EMR t in function of W, NemHi 
and, hence, we can calculate the ratio of the rates as a func- 
tion of the inclination l and spin s: 



Kerr Schw , — t — / \ 

flEMRI = CSEMRI X W 6 - 2 '' ((,, S) 



x VV 12 - 4 ^ (i, s) . 



(34) 
(35) 



For instance, for a spin of s = 0.999 and an inclination of i — 

0. 4 rad, we estimate that W ~ 0.26 and, thus, Nemri ~ 3 °. 

1. e. we boost the event rates by a factor of 30 in comparison 
to a non-rotating MBH. 



B EMRI — K 1 R 

K- 



1 M. 



In A N.m. / 

[i6.977y (3 - 7) (1 + 7 ) 3/2 ] . (30) 

We can then derive the event rate for the Schwarzschild case, 
based on (171 



A& * 0.235 4 ( 3 -?!^ 1+ ^ /2 



In A 6-2^ 



In 



9-4 7 

N.m. \ 
M. 



GM. 



(31) 



The last equation is based on the assumption that equa- 
tion (27l holds, and the last equality, (1 — e)a — 8GM,/c 2 , 
is the "effective" value of the periapsis for the last parabolic 
stable orbit (LSO from now onwards) for a Schwarzschild 
MBH. I.e. if the star is on an orbit that in Newtonian dy- 
namics leads to a periapsis smaller than this value, the star 
disappears if we take into account relativistic dynamics. 



6.2 The Kerr case 

In the Kerr case we simply have to recalculate where this 
LSO is by taking into account the value of the spin of the 



7 NET EFFECT OF RESONANT 
RELAXATION ON EMRI RATES 

7.1 The role of vectorial resonant relaxation 

To understand the impact of the previous calculation on 
event rates for EMRIs we have to elaborate on prograde and 
retrograde orbits. We have seen that retrograde EMRI orbits 
see the central MBH as an effectively "larger" MBH; i.e. it 
is easier to plunge through the horizon. Contrary, prograde 
EMRI orbits "see" the effective size of the MBH shrink and, 
thus, have a harder time in hitting the central MBH. It is 
therefore important to assess the orientation of orbits in the 
regime of interest. It takes on average (vectorial) resonant re- 
laxation (RR) a time t RRj v to rotate coherently the orbital 
plane of an orbit by an angle it/ 2 (Hopman & Alexander 



2 



2006 [ ). To change a prograde (retrograde) orbit to a retro- 
grade (prograde) orbit, it takes four times longer: The 7r/2 
rotation is the maximum that can be obtained over the self- 
quenching time; the rest to get up to a full n rotation is 



done non-coherently over 4 coherence times (see Bregman & 



Alexander||2011 for a discussion of the numeric prefactors). 



It should be noted that vector RR is invariant under 
precession (see e.g. |Hopman fc Alexan der 2006). We must 
note also that the change in the inclination of the orbit with 
respect to the spin axis due to GW emission is relatively 



© 2012 RAS, MNRAS 000, [TV?? 



Spin matters: EMRI event rates 9 



Spin (s) 


Inclination (t, rad) 


Vv\i. s ) 


0.900 


0.100 


0.429452 


0.900 


0.400 


0.448093 


0.900 


0.700 


0.499450 


n onn 
u.yuu 


1 onn 
l.UUU 


U.Oyoz / o 


n onn 
u.yuu 


1 Qnn 
l.oUU 


n T'QOQQO 
U. / 013001) 


n onn 
u.yuu 


1.0 ( u 


U.ooou ( y 


0.900 


-0.100 


1.415955 


0.900 


-0.400 


1.377239 


0.900 


-0.700 


1.295011 


n onn 
u.yuu 


-l.UUU 


1.1 to 1 UU 


0.950 


0.100 


0.370036 


0.950 


0.400 


0.386009 


0.950 


0.700 


0.436921 


n ot^n 

u.you 


1 onn 
l.UUU 


U.04oo0z 


u.you 


i Qnn 
l.oUU 


U. / UoZO 1 


n oc;n 

u.you 


1.0 (U 


n q^vqoo 
U.oO ( ozU 


0.950 


-0.400 


1.396449 


0.950 


-0.700 


1.309052 


0.950 


-1.000 


1.181942 


U.yDU 


1 Qnn 
-l.oUU 


i.U^4oOD 


0.990 


0.100 


0.297301 


0.990 


0.400 


0.306924 


0.990 


0.700 


0.354716 


0.990 


1.000 


0.494738 


n oon 
u.yyu 


1 Qnn 
l.oUU 


U.O ( y4Do 


n oon 
u.yyu 


1.0 I u 


U.0OZ0ZI 


0.990 


-0.100 


1.454732 


0.990 


-0.400 


1.411720 


0.990 


-0.700 


1.320145 


0.990 


-1.000 


1.186631 


u.yyu 


i Qnn 
-l.oUU 


1 nonQ 1 yi 
!.UzUol4 


0.999 


0.100 


0.260205 


0.999 


0.400 


0.264063 


0.999 


0.700 


0.310302 


0.999 


1.000 


0.479038 


0.999 


1.300 


0.672349 


0.999 


1.570 


0.849364 


0.999 


-0.100 


1.458589 


0.999 


-0.400 


1.415145 


0.999 


-0.700 


1.322624 


0.999 


-1.000 


1.187655 


0.999 


-1.300 


1.019828 



Table 2. Values for W. 



rather small (see [Hughes (20 00] |200Tj >), so small that fre- 
quently it has been assumed to be constant, which provides 
an extra equation for the evolution of the Carter constant 
in the inspiral process, making things significantly simpler. 

The dependence of the transverse RR torque (i.e. 
direction-changing torque) on the eccentricity has been mea- 
sured from simulations by Giirkan & Hopman (20071. In 



their work, the authors derive that it grows quadratically 
by a factor 3 in total from to 1. 

The radius of the sphere of influence is 



1 GM. , 1 

finfl = — — —T~ ~ 1 P C 1 



1+7 < 



M. 



1+7 V 10 6 M Q 



60 km/s 



(Jo 



(36) 

for a given exponent 7. Within this radius the relaxation 
time is 



/ x N-3/2 lnA<r 3 (r) 

t r (r oc (l + 7 ) 3/2 --— 

G 2 {m)m co n(r) 

~ 2 x 10 8 yr(l + 7r 3/2 ( 



100 kms" 1 



10 M, 



10 6 M Q pc~ 3 
(m)n 



(37) 



In this equation we follow the usual notation: o(r) is the 
local velocity dispersion; for r < r infl it is approximately the 
Keplerian orbital speed \/GM.i — 1 ; n(r) is the local density 
of stars, (m) is the average stellar mass, m co is the indi- 
vidual mass of the compact object, which we assume to be 
all SBHs, and take a mass of m co = 10 M Q for all of them. 



In the vicinity of a MBH, (r < r infl ), A w M./m* ( Bah- 
call fc Wolf|1976||Lightman fc Shapiro|1977) , and typically 
InA ~ 11. 

Relaxation redistributes orbital energy amongst stellar- 
mass objects until SBHs form a power-law density cusp, 
n(r) oc r ' with 7 ~ 1.75 around the MBH, while less mas- 
sive species arrange themselves into a shallower profile, with 
ct — 1.4 — 1.5 as we have mentioned earlier, although recent 
studies have found a general solution to the problem of mass 
segregation around MBH in galactic nuclei, with a more ef- 
ficient diffusion for the heavy stars, reaching a 7 ~ —2 in 



the "strong 


; mass-segregation" regime (Alexander & Hop- 


man||2009| 


Preto & Amaro-Seoanc 2010| |Amaro-Seoane & 


Preto|2011 


1. 



Since cr(r)~ = GM,/r and we take that (mi) = 0.7 A/p and, 
as mentioned, mco — 10 -^o 1 we have all information to 

, peri(^) 



derive t T pcli (r) := (1 — e)t T from Equation ( |37[ ) and (18 1. 

As regards the explicit expression for the characteristic 
timescale for vectorial resonant relaxation, from |Hopman &1 
Alexander ( 2006 1 we have that 



P(a 



M, (3 v (ef 



(38) 



where we have taken into account the corrections for high 



values of e as given in Giirkan & Hopman (20071, f3 v {e) = 
0.28 (e 2 + 0.5) and P(a) = a''7(GM,) 1/2 . This allows us to 
follow the dependence with the radius (and eccentricity) of 



the ratio t r 



r . pen 



If we now equate the timescales of interest, the grav- 
itational radiation driven time t GW , defined as in the ap- 



proximation of Peters ( 1964 ) , to the two-body relaxation 



time at periapsis, t l poli , we obtain the short-dashed curve 
of figure [3] on the left of this line, the contribution of GW 
radiation to orbital evolution dominates over two-body re- 
laxation. In the absence of resonant relaxation, if a SBH 
crosses this line from the right (lower eccentricities), it will 
become an EMRI, provided, of course, that it is still on a 
stable orbit, i.e. above the separatrix corresponding to its 
orbital orientation. For a Schwarzschild MBH, all separatri- 
ces are the same and there is a unique critical point (PS). A 
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SBH with a semi-major axis larger that the value at PS will 
experience a direct plunge if relaxation brings its eccentric- 
ity to a high value because it will cross the separatrix (and 
be swallowed in less than an orbital time) before it has a 
chance to enter into the GW-dominated regime. Conversely, 
objects with smaller semi-major axis values are much more 
likely to end up as EMRIs rather than plunges. 

For a fast spinning SMBH, the separatrix for prograde 
orbits is shifted to significantly lower a values, with a corre- 
sponding higher value of the critical semi-major axis, corre- 
sponding to the point PP in figure [3] As we have explained 
above, it is this effect which can lead to a significant increase 
in the EMRI rate, combined with the fact that the critical 
point for retrograde orbits (PR) is much less affected and 
that an isotropic orbit distribution is expected, thanks to 
relaxational processes. However this increase in EMRI rate 
would can be thwarted by vector RR if this process can 
change the orbital orientation of a SBH after it has crossed 
the "iGw — ^r.pcri" li ne an d before it has completed its GW- 
driven inspiral, i.e. on a timescale shorter than t Gw . Indeed, 
if the orbit becomes significantly less prograde as the the in- 
spiral takes place, due to RR, the separatrix moves up and 
the SBH might suddenly find itself on a plunge orbit. 

To check for this possibility, we also plot, in figure [3j 
a long-dashed line corresponding to the condition t GW = 
( RR „ with t GW < tun v on the left of this line. SBHs that 
cross the u t GW — f r ,pcri" line while on the left side of the 
"t GW = i RR v " line keep their orbital orientation during 
their inspiral and complete it without abrupt plunge. One 
can see that, for our choice of parameters, this is the case 
for all prograde orbits. On the other hand, retrograde orbits 
can cross the "tew = ^r,pcri" une while RR is still effective 
enough to change their orientation during inspiral. However, 
the effect of RR on retrograde orbits cannot reduce signifi- 
cantly the total EMRI rate and may even increase it slightly 
because (1) these orbits contribute less that the prograde 
ones (and more to the plunge rate) and (2) statistically, RR 
is more likely to make the orbit become less retrograde which 
pushes down the separatrix. 

Finally, we also note that for the other proposed mecha- 
nism to produce EMRIs, the tidal separation of a binary con- 
taining a compact object ( Miller et al.|2005 Amaro-Seoane 



et al.|2007 l, the captured objects typically have much lower 



eccentricities and smaller semi-major axis. Therefore, they 



cross \ w = t T 



line and start their inspiral, with or- 



bital parameters well above the uppermost separatrix (for 
retrograde orbits). As the GW-driven trajectory in the e — a 
plane is basically parallel to the separatrices, there is no 
danger of a premature plunge, even though RR has ample 
time to flip the orbital orientation during inspiral. 



EMRIs, as a number of different studies by different authors 
using different methods find. 

Up to now spin effects of the central MBH have been 
always ignored. Hence, the question arises, whether a plunge 
is really a plunge when the central MBH is spinning. This 
consideration has been so far always ignored. 

So as to estimate EMRI event rates, one needs to know 
whether the orbital configuration of the compact object is 
stable or not, because this is the kernel of the difference be- 
tween an EMRI and a plunge. In this paper we take into 
account the fact that the spin makes the LSO to be much 
closer to the horizon in the case of prograde orbits but it 
pushes it away for retrograde orbits. Since the modifications 
introduced by the spin are not symmetric with respect to 
the non-spinning case, and they are more dramatic for pro- 
grade orbits, we prove that the inclusion of spin increases 
the number of EMRI events by a significant factor. The ex- 
act factor of this enhancement depends on the spin, but the 
effect is already quite important for spins around s ~ 0.7. 

We also prove that these fake plunges, "our" EMRIs, 
do spend enough cycles inside the band of eLISA to be de- 
tectable, i.e. they are to be envisaged as typical EMRIs. We 
note here that whilst it is true that EMRIs very near the 
new separatrix shifted by the spin effect will probably con- 
tribute not enough cycles to be detected, it is equally true 
for the old separatrix (Schwarzschild, without spin). In this 
sense, we find that the spin increases generically the number 
of cycles inside the band for prograde EMRIs in such a way 
that EMRIs very near to the non-spin separatrix, which con- 
tributed few cycles, become detectable EMRIs. In summary, 
spin increases the area, in configuration space of detectable 
EMRIs. We predict thus that EMRIs will be highly domi- 
nated by prograde orbits. 

Moreover, because spin allows for stable orbits very near 
the horizon in the prograde case, the contribution of each 
cycle to the SNR is significantly bigger than each cycle of 
an EMRI around a non-spinning MBH. 

We then show that vectorial resonant relaxation will 
not be efficient enough to change prograde orbits into ret- 
rogrades once GW evolution dominates (which would make 
the EMRIs plunge instantaneously, as they would be in a 
non-allowed region of phase space). 

These new kind of EMRIs we describe here, high- 
eccentric EMRIs, are produced by two-body relaxation and, 
as such, they are ignorant of the Schwarschild barrier. While 
low-eccentricity EMRIs run into the problem of having to 
find a way to travers this barrier, our "plunge-EMRIs" do 
not. We predict that EMRI rates will be dominated by high- 
eccentricity binaries, with the proviso that the central MBH 
is Kerr. 



8 CONCLUSIONS 

In this article we have addressed the problem of direct 
plunges and MBHs. If this is a Schwarzschild MBH, the com- 
pact object will plunge through the horizon and will hence 
not contribute to the mapping of space and time around the 
MBH, contrary to an EMRI, which describes thousands of 
cycles before it merges with the central MBH. On the other 
hand, the event rate of plunges is much larger than that of 
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